Effective mode volumes and Purcell factors for leaky optical cavities 
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We show that for optical cavities with any finite dissipation, the term "cavity mode" should 
be understood as a solution to the Helmholtz equation with outgoing wave boundary conditions. 
This choice of boundary condition renders the problem non-Hermitian, and we demonstrate that 
the common definition of an effective mode volume is ambiguous and not applicable. Instead, we 
propose an alternative effective mode volume which can be easily evaluated based on the mode 
calculation methods typically applied in the literature. This corrected mode volume is directly 
applicable to a much wider range of physical systems, allowing one to compute the Purcell effect 
and other interesting optical phenomena in a rigorous and unambiguous way. 



Optical microcavities are inherently dissipative and are 
typically characterized by a quality factor, or Q-value, 
describing the relative energy loss per cycle as well as 
an effective mode volume, T4ff, which gives a measure 
of the spatial confinement of the electromagnetic field in 
the cavity. Cavities with high Q-values and small mode 
volumes provide enhanced light-matter interaction and 
are of fundamental as well as technological interest [l|, y] . 
Effective mode volumes are ubiquitous in physics and 
connect to a wide range of phenomena, including sens- 
ing [3|], switching [J], cavity quantum electrodynamics 
(QED) 0, circuit-QED 0, and optomechanics Q- As a 
striking example of the use of mode volumes, an emitter 
in an optical cavity will experience a medium-enhanced 
radiation rate relative to that in a homogeneous medium 
given by the so-called Purcell factor [8[ 
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where Ac is the free space wavelength, and ric is the ma- 
terial refractive index at the field antinodc r^. Purcell 
factors are widely used in quantum optics as a figure of 
merit for single photon sources |9|. 

Mode volumes are often attributed to the physically 
appealing idea of a single cavity mode. However, in spite 
of the fact that cavity modes are widely used in the lit- 
erature, there seems to be a disturbing lack of a precise 
definition, and their mathematical properties therefore 
remain unspecified. The lack of a definition is evidenced 
in part by the diverse nomenclature at use ( "resonance" , 
"leaky mode" or "quasi mode"), suggesting that the dis- 
sipative nature of cavity modes somehow makes them 
different from other modes, but an explicit distinction is 
rarely made. It appears that cavity modes are widely be- 
lieved to share the properties of Hermitian eigenvectors, 
although a mode with a finite lifetime is incompatible 
with the solution space of Hermitian eigenvalue prob- 
lems. Strictly, only in the limit of large (infinite) Q will 
the modes in optical cavities appear as the solutions to 
Hermitian eigenvalue problems. For finite Q, the lack of 
hermiticity effectively renders expressions such as Eq. ([T]) 



ambiguous, since the volume VeS cannot be inferred from 
the usual inner product of Hermitian systems. In particu- 
lar, if er(r) describes the relative permittivity distribution 
and fc (r) is the cavity mode, then a direct application of 
the common normal mode prescription 
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cannot be expected to provide the correct mode volume. 
In many practical calculations, the leaky cavity mode 
is found with the finite-difference time-domain (FDTD) 
method by launching a short pulse and monitoring the 
resonant field that leaks from the cavity at a rate set by 
the Q-value [10[ . Figure [1] shows a sketch of an example 
cavity along with mode profiles calculated with FDTD 
111 . For this cavity mode the integral in Eq. ([2]) diverges 
as a function of the integration volume V , and indeed 
this is formally the case for all cavities with a finite Q. 
For very high-Q cavities, however, the divergence is slow 
and may not be discernable in practice due to numerical 
accuracy, but the formal divergence still renders Eq. ^ 
questionable. 




FIG. 1. Sketch of a photonic crystal made from a triangular 
lattice of air holes (lattice constant a) in a membrane of high 
refractive index. A defect cavity is formed by the omission of 
a single hole. Right: Absolute value of the t/-polarized cavity 
mode in the planes z = Q (top) and y = (bottom). 

In this Letter we argue that the term "cavity mode" 
should be understood as a solution to the Helmholtz 



equation with outgoing wave boundary conditions. This 
definition renders the cavity modes identical to the quasi- 
normal modes of Lee et al. [l^ which have complex res- 
onance frequencies (with a negative imaginary part, as 
expected for dissipative modes) and exhibit an inherent 
exponential divergence at large distances. We illustrate 
directly how this definition complies with typical calcula- 
tions of cavity modes using FDTD and we elucidate how 
the cavity modes from FDTD calculations also show an 
exponential divergence. Quasinormal modes have prop- 
erties that are different from the solutions to Hermitian 
eigenvalue problems and therefore important results such 
as orthogonality and completeness of the solutions can- 
not be taken for granted. Nevertheless, an inner product 
and a corresponding orthogonality relation can be de- 
fined, and the quasinormal modes can be used as a basis 
for expansion of the electromagnetic Green's tensor in 
certain regions [12|. This enables a precise and unam- 
biguous description of light-matter interaction in general 
leaky optical cavities, including an expression for the ef- 
fective mode volume and thus the Purcell factor. 

It is illustrative to start by highlighting the differ- 
ences between two types of modes that can be associ- 
ated with fields in optical cavities. The electric field in 
general non-magnetic materials satisfy the wave equa- 
tion with time-harmonic solutions of the form E(r,i) = 
E(r, oj) cxp{— iwi}. The position-dependent electric field 
E(r) solves the vector Hclmholtz equation 



normalized as llj 



V X V X E(r,w) - A:^e,.(r)E(r,w) = 0, 



(3) 



where fco = lo / c. Together with a suitable set of bound- 
ary conditions, Eq. ([3|) provides a generalized eigenvalue 
equation. We will use the term normal mode to denote 
a solution to Eq. ([3]) with any set of boundary condi- 
tions that renders the problem Hermitian. In this case 
we denote the vector eigenfunctions and corresponding 
real cigcnfrcquencies as f^(r) and w^^, respectively. The 
normal modes are normalized as 



(fJfA)= /e.(r)f;(r).f;,(r)dr 
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where the integral is over the volume defined by the 
boundary conditions. In many applications, the limit 
y — )• cxD is taken, in which case the spectrum of eigenval- 
ues becomes continuous. We will use the term quasinor- 
m,al m,odes for solutions to Eq. ([3]) with outgoing wave 
boundary conditions (the Sommerfeld radiation condi- 
tion [l3|)- This choice of boundary condition renders 
the eigenvalue problem non-Hermitian with a discrete 
spectrum, and we denote the vector eigenfunctions with 
a tilde as f/i(r). The corresponding eigenfrequencies, 
' iwj,, are in general complex with wj^ < 0, and 
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it follows from Eq. ^ that, contrary to the Hermitian 
case, f^(r) and f*(r) are not eigenvectors corresponding 
to the same eigenvalue. The quasinormal modes may be 



((fJfA»= lim / er(r)f^(r).f;,(r)dr 
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where dV denotes the border of the volume V. The limit 
T^ — >■ oo is calculated by increasing the volume to obtain 
convergence. For the systems that we investigate in this 
Letter, the convergence is remarkably fast. For very low- 
Q cavities, however, the convergence is nontrivial due 
to the exponential divergence of the quasinormal modes 
which may cause the inner product to oscillate around 
the proper value as a function of calculation domain size. 
In addition to the modes of the cavity it is convenient 
to introduce the electromagnetic Green's tensor through 

VxVxG(r,r',cj)-fc2e,(r)G(r,r',a;) = K(r-r'), (6) 

subject to the Sommerfeld radiation condition. The 
Green's tensor provides the proper framework for cal- 
culating light emission and scattering in general dielec- 
tric structures. In general, the decay rate r„(r,w) of a 
dipole emitter with orientation e^ may be enhanced or 
suppressed as compared to the rate Fb in a homogeneous 
medium. The relative rate may be expressed as 

Ta{r,uj) _ lm{eaG{r,r,uj)ea} 



Tb{uj) Im{eaGB(r,r,a;)eQ}' 
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where Gb (r, r', oj) is the Green's tensor in a homogeneous 
background medium with ei(r) = eb |14| . In certain re- 
gions, such as inside the scattering region, the transverse 
part of the Green's tensor may be expanded through [12 1 
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The implicit assumption behind the notion of a cavity 
mode is that one term ^ = c dominates the expansion 
in Eq. (|8]) and hence that the Green's tensor can be ap- 
proximated by this term only. With this assumption, 
and noting that Im{G(r, r,a;)} = Im{G'^(r, r,a;)}, one 
can use Eqs. ([7]) and ^ with u! = uic = ^f to recover 
Eq. ([1} with a corrected effective mode volume given as 



K 



Q 



1 \VQ? 



eff 



«Q 



((fclfc)) 

fc^(rc) 



(9) 



where 



«Q 



Uq + ivQ is complex in general. This pre- 
scription provides a direct and unambiguous way of cal- 
culating the effective mode volume. 

The quasinormal modes can be calculated analytically 
for sufficiently simple structures, but for general struc- 
tures the outgoing wave boundary conditions are not im- 
mediately compatible with standard numerical solution 
methods. One option is to rewrite Eq. ([3]) as 

V X V X E(r,u;) - A:^eBE(r,a;) = /c^A£(r)E(r, w), (10) 



where Ae — er(r) — ee, and calculate the quasinormal 
modes from a Fredholm type integral equation, 

E(r,c.) = (^)Y G^(r,r',c.)A£(r')E(r',c.)dr', (11) 

which manifestly respects the outgoing wave boundary 
conditions. Another option is to use FDTD to calculate 
the quasinormal mode as the resonant field that is ex- 
cited by an initial short input pulse. We have used both 
Eq. PT|) and FDTD to calculate the quasinormal modes 
in different example cavities in two and three dimensions. 
For our example below, we solve Eq. (|TT|) using the ex- 
pansion technique of Ref. |15| with an additional iteration 
of fco to make the solution self-consistent. In addition, we 
perform FDTD calculations using perfectly matched lay- 
ers (PMLs) [16[ to enforce the outgoing wave bounday 
conditions. 

We first consider a 2D finite-sized hexagonal crystal- 
lite of high-index rods in air with a single missing rod in 
the center. The rods have relative permittivity £1=11.4 
and radius i?=0.15 a, where a is the lattice constant. We 
focus on transverse magnetic (TM) polarization in which 
the electric field is in the direction of the rods. In the 
limit of infinite s ize, the photonic crystal exhibits a pho- 
tonic band gap [l7|, and the Q of the cavity therefore 
depends on the size of the structure which we may char- 
acterize by the number of rod layers, N. For iV=l and 
N=2, Fig. [2] shows the supported cavity modes with fre- 
quencies a)ca/27rc=0. 4259— 0.0135i {Q = —Cjyi/ICji ~ 16) 
and Wca/27rc = 0.4218 -0.0013z {Q « 163). In the case of 
A^=3 (not shown), the structure supports a cavity mode 
with frequency 0j^a/2TTc = 0.4216 - O.OOOli [Q k. 1576). 
As expected, the quasinormal modes are concentrated 
in the center of the cavity and seem to fall off with in- 
creasing distance to the crystallite. At large distances, 
however, the quasinormal modes (by definition) behave 
as outgoing waves of the form f(r) ex exp(zA:o?')/-\/?' (2D) 
and f (r) oc eyi^{ikQr)/r (3D), and since A:o = fcn + iki 
with fci < 0, they diverge exponentially as r ^ oo. For 
the case of 7V=1, the top panel in Fig. [2] illustrates di- 
rectly that the solutions to Eq. ([TT]) are identical to those 
obtained from FDTD. In particular, both solution meth- 
ods pick up the divergence in the field at large distances. 
Figure [3] shows, as a function of the size of the calcu- 
lation domain, the corrected effective mode volume in 
Eq. © along with the common definition in Eq. ^ with 
fc(r) = fc(r). From the figure it is clear that whereas 
V^ converges quickly to the limiting values, V^ seems 
to increase with the size of the domain. The linear diver- 
gence in V^g with the size of the normalization domain 
was also noted in Ref. [18| and derives from the fact that 
the field does not go to zero at positions outside the crys- 
tallite, cf. Fig. [21 At much larger V, the field, and hence 
V^, diverges exponentially. For increasing Q, the linear 
divergence with domain size becomes less and less pro- 
nounced, suggesting that the two formalisms converge in 
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FIG. 2. (a): Absolute value along the x-sxis of the quasinor- 
mal mode in the 2D crystallite for the case of A''=l. Blue solid 
line shows the solution to Eq. (|11|) . and black circles show the 
calculation using FDTD. Inset shows long distance behavior 
on a logarithmic scale, (b): Absolute value along the x-axis 
of the quasinormal mode for the case of A^ = 2 with the inset 
showing the distribution in the sy-plane. Grey shaded areas 
indicate the high-index rods. 



the limit of infinite Q as expected. 

Next, for a practical 3D example we consider a pho- 
tonic crystal membrane (ei.=12) of thickness h = 0.5 a 
and hole radius r = G.275a. A single air hole is omit- 
ted to create a cavity, and Fig. [1] shows the supported 
cavity mode with frequency (ija/27rc=:0.2904 — 0.0004i 
{Q ~ 362). From the results in two dimensions we know 
that the quasinormal modes can be directly calculated 
using FDTD with PMLs. The field in Fig. [U was cal- 
culated in the same way, and we therefore argue that 
it is indeed a quasinormal mode and that we should use 
Eq. ^ rather than Eq. ^ to calculate the effective mode 
volume. Fig. [4] shows both V^ and V^ as a function of 
calculation domain size. At the quasinormal mode fre- 
quency, the photonic band gap prevents in-plane propa- 
gation, and therefore the only way for the field to leak out 
of the cavity is in the z-direction. This means that both 
V^^ and V^ converge quickly as a function of width and 
depth of the calculation domain and we focus only on the 
variation in the effective mode volumes with the height 
of the calculation domain. As in two dimensions, the 
data shows a fast convergence of V^, while V^ clearly 
diverges, confirming that Eq. ^ is not applicable. 
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FIG. 3. EfTective mode volumes V^ (thick lines) and V^ 
(thin lines) for N — 1 (red dash-dotted), N — 2 (green 
dashed) and A'^ = 3 (blue solid) as a function of radius R 
of the calculation domain. Circles indicate reference mode 
volumes V^g^ from independent Green's tensor calculations 
[la ], and grey dashed areas show the rod cross sections along 
the a;-axis. 
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FIG. 4. Effective mode volume V^ (red dashed) and V^ 
(blue solid) for the cavity in Fig. [l] as a function of height 
of the calculation domain. Circles indicate reference mode 
volumes Vj^ derived from independent Green's tensor cal- 
culations [10] with estimated error bars at different domain 
heights. Gray dashed area shows the extend of the membrane. 



Finally, we compare the calculated mode volum.es to in- 
dependent calculations using the Green's tensor [1^ [l5 ] . 
Substituting Fp = rc(r,a;c)/rB(a;c) in the expression for 
the Purcell factor [19| defines an efi^ective mode volume 
V^*^*. For each of the cavities, V^g' is indicated with a 
circle in Fig. [3l The maximum estimated absolute error 
in these calculations is less than 0.0003. The observable 
discrepancies for A^=l and N=2 stem from the single 
mode approximation and indicates the limited validity of 
the Purcell factor. In Fig. [JJ V^f^^ was calculated with 
FDTD as the response to an input dipole source at three 
different domain sizes and with estimated error bars as 
indicated. These independent calculations confirm that 
Eq. ((9|) not only is unambiguous, but also leads to the 
correct value within the single mode approximation. 

In conclusion, we have shown that the term "cavity 



mode" should be understood as a so-called quasinor- 
mal mode, defined as a solution to the Helmholtz equa- 
tion with outgoing wave boundary conditions. This can 
have profound consequences, since this choice of bound- 
ary conditions renders the differential equation problem 
non-Hermitian so that common results from Hermitian 
eigenvalue analysis do not apply. In particular, the quasi- 
normal modes have complex frequencies and exhibit an 
inherent divergence at long distances which makes the 
calculation of an effective mode volume nontrivial. In- 
troducing an inner product that carefully accounts for 
the long distance behavior, it is possible to normalize the 
quasinormal modes and to define an effective mode vol- 
ume in a direct and unambiguous way. In practical cal- 
culations, this corrected mode volume can be obtained 
in a straightforward way using exactly the same cavity 
modes that are typically computed for use in mode vol- 
ume calculations. 
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